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Tuning the magnetism of ordered and disordered strongly-correlated 

electron nanoclusters 

Nicholas Kioussis, Yan Luo, and Claudio Verdozzjj 
Department of Physics, California State University Northridge, California 91330-8268 

Recently, there has been a resurgence of intense experimental and theoretical interest on the 
Kondo physics of nanoscopic and mesoscopic systems due to the possibility of making experiments 
in extremely small samples. We have carried out exact diagonalization calculations to study the 
effect of energy spacing A in the conduction band states, hybridization, number of electrons, and 
disorder on the ground-state and thermal properties of strongly-correlated electron nanoclusters. For 
the ordered systems, the calculations reveal for the first time that A tunes the interplay between 
the local Kondo and non local RKKY interactions, giving rise to a "Doniach phase diagram" for the 
nanocluster with regions of prevailing Kondo or RKKY correlations. The interplay of A and disorder 
gives rise to a A versus concentration T = phase diagram very rich in structure. The parity of the 
total number of electrons alters the competition between the Kondo and RKKY correlations. The 
local Kondo temperatures, Tk, and RKKY interactions depend strongly on the local environment 
and are overall enhanced by disorder, in contrast to the hypothesis of "Kondo disorder" single- 
impurity models. This interplay may be relevant to experimental realizations of small rings or 
quantum dots with tunable magnetic properties. 



I. INTRODUCTION 



Magnetic impurities in non-magnetic hosts have been 
one of the central subjects in the physics of strongly cor- 
related systems for the past four decadesi. Such endur- 
ing, ongoing research effort is motivated by a constant 
shift and increase of scientific interest over the years, from 
dilute^ to concentrated impurities'^ , from periodic^ to dis- 
ordered samples-'^, and from macroscopic^ to nanoscale 
phenomena™. Macroscopic strongly correlated electron 
systems at low temperatures and as a function of mag- 
netic field, pressure, or alloying show a wide range of 
interesting phenomena, such as non-Fermi-liquid behav- 
ior, antiferromagnetism, ferromagnetism, enhanced para- 
magnetism, Kondo insulators, quantum criticality or 
superconductivitjiiii. These phenomena are believed to 
arise through the interplay of the Kondo effect, the elec- 
tronic structure and intersite correlations. In the sim- 
plest single-impurity case, the Kondo problem describes 
the antiferromagnetic interaction, J, between the impu- 
rity spin and the free electron spins giving rise to an 
anomalous scattering at the Fermi energy, leading to a 
large impurity contribution to the resistivitj*^ . The low- 
energy transport and the thermodynamic properties scale 
with a single characteristic energy, the Kondo tempera- 
ture, Tk oc exp{—l/p{EF)J), where p{Ep) is the den- 
sity of states of the conduction electrons at the Fermi 
energjii. At T >> Tk, the impurity spin is essentially 
free and the problem can be treated perturbatively. At 
T « Tk , the impurity spin is screened forming a singlet 
complex with the conduction electrons, giving rise to a 
local Fermi liquid state. 

For dense Kondo or heavy fermion compounds con- 
taining a periodic array of magnetic ions interacting with 
the sea of conduction electrons, the physics is determined 
from the competition between the non local Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interactions and the local 



Kondo interactions^. The RKKY interaction is an in- 
direct magnetic interaction between localized moments 
mediated by the polarized conduction electrons, with an 
energy scale of order Jrkky oc J'^p{Ef), which pro- 
motes long- or short-range magnetic ordering. On the 
other hand, the Kondo effect favors the formation of 
singlets resulting in a non-magnetic ground state. In 
the high temperature regime the localized moments and 
the conduction electrons retain their identities and inter- 
act weakly with each other. At low-temperatures, the 
moments order magnetically if the RKKY interaction is 
much larger than the Kondo energy, while in the reverse 
case, the system forms a heavy Fermi liquid of quasipar- 
ticles which are composites of local moment spins bound 
to the conduction electrons^iS. Thus, the overall physics 
can be described by the well-known "Doniach phase dia- 
gram", originally derived for the simple Kondo necklace 
modeliS.. The description of the low-temperature state, 
when both the RKKY and the Kondo interactions are 
of comparable magnitude, is an intriguing question that 
remains poorly understood and is the subject of active 
research^. 

The interplay of disorder and strong correlations has 
been a subject of intensive and sustained research, in 
view of the non- Fermi-liquid (NFL) behavior in the vicin- 
ity of a quantum critical poinliii. Various disorder-driven 
models have been proposed to explain the experimen- 
tally observed^ NFL behavior at low temperatures^iSiLiS. 
The phenomenological "Kondo disorder" approaches^ii^, 
based on single-impurity models, assume a distribution 
of Kondo temperatures caused by a distribution of either 
f — c orbital hybridization or of impurity energy levels. 
These models rely on the presence of certain sites with 
very low Tk spins leading to a NFL behavior at low T. 
An open issue in such single-site Kondo approaches, is 
whether the inclusion of RKKY interactions would renor- 
malize and eliminate the Iow-Tr- spin3^ii^ii^*i&. An alter- 
native view is the formation of large but finite magnetic 



clusters (Griffith phases) within the disordered phase 
through the competition between the RKKY and Kondo 
interactionsfiiii. 

On the other hand, the possibihty of making exper- 
iments in extremely small samples has lead to a resur- 
gence of both experimental and theoretical interest of 
the physics of the interaction of magnetic impurities in 
nanoscopic and mesoscopic non-magnetic metallic sys- 
tems. A few examples include quantum dotsi^, quantum 
boxegiS. and quantum corrals^. Recent scanning tun- 
neling microscope(STM) experiments^^ studied the in- 
teraction of magnetic impurities with the electrons of a 
single-walled nanotube confined in one dimension. Inter- 
estingly, in addition to the bulk Kondo resonance new sub 
peaks were found in shortened carbon nanotubes, sepa- 
rated by about the average energy spacing. A, in the nan- 
otube. The relevance of small strongly correlated systems 
to quantum computation requires understanding how the 
infinite-size properties become modified at the nanoscale, 
due to the finite energy spacing A in the conduction 
band which depends on the size of the particl924i2422iS4S4. 
For such small systems, controlling Tr- upon varying 
A is acquiring increasing importance since it allows to 
tune the cluster magnetic behavior and to encode quan- 
tum information. While the effect of A on the single- 
impurity Anderson or Kondo model has received con- 
siderable theoreticalSii2iS2iSiSi and experimentalSli atten- 
tion recently, its role on dense impurity clusters remains 
an unexplored area thus far. The low-temperature be- 
havior of a nanosized heavy-electron system was recently 
studied within the mean-field approximation^^. A cen- 
tral question is what is the effect of A on the interplay 
between the Kondo effect and the RKKY interaction 

In this work we present exact diagonalization calcula- 
tions for d- or /-electron nanoclusters to study the effects 
of energy spacing, parity of number of electrons, and hy- 
bridization on the interplay between Kondo and RKKY 
interactions in both ordered and disordered strongly cor- 
related electron nanoclusters. While the properties of the 
system depend on their geometry and sizc^^, the present 
calculations treat exactly the Kondo and RKKY interac- 
tions, the disorder averages, and they provide a distribu- 
tion of local Tk's renormalized by the intersite f-f interac- 
tions. Our results show that: i) tuning A and the parity 
of the total number of electrons can drive the nanocluster 
from the Kondo to the RKKY regime, i.e. a zero- temper- 
ature energy spacing versus hybridization phase diagram; 
ii) the temperature versus hybridization "Doniach" phase 
diagram for nanoclusters depends on the energy spac- 
ing ; iii) changing the total number of electrons from 
even to odd results in an enhancement (suppression) of 
the local Kondo (RKKY) spin correlation functions; iv) 
the A versus alloy concentration T — phase diagram 
exhibits regions with prevailing Kondo or RKKY corre- 
lations alternating with domains of ferromagnetic (FM) 
order; and v) the local Tr-'s and the nearest-neighbor 
(n.n) RKKY interactions depend strongly on the local 
environment and are overall enhanced by disorder. The 



disorder-induced enhancement of Tk in the clusters is in 
contrast to the hypothesis of "Kondo disorder" models 
for extended systems. 

The rest of the paper is organized as follows. In Sec. 
II, we describe the model for both the periodic and disor- 
dered clusters. In Sees. IIIA and IIIB we present results 
for the ground-state and thermal properties of the or- 
dered and disordered nanoclusters, respectively. Section 
IV contains concluding remarks. 



II. METHODOLOGY 

The one dimensional Anderson lattice model is 

H = -i ^(cLc,+i, + H.c) + EfJ2 4a 

i(7 i<7 

+uY.<Al +Vj2iflc,. + H.C.). (1) 

i ia 

Here, t is the nearest-neighbor hopping matrix element 
for the conduction electrons, cf^{ci^a) and f^_„{fi,a) cre- 
ate (annihilate) Wannier electrons in c- and /- like orbital 
on site i with spin a, respectively; Ef is the energy level 
of the bare localized orbital, V is the on-site hybridiza- 
tion matrix element between the local f orbital and the 
conduction band and U is the on-site Coulomb repul- 
sion of the f electrons. We use a simple nearest-neighbor 
tight-binding model for the conduction band dispersion, 
ek = -2tcosk. We consider the half-filled {Nei = 2N) 
symmetric {Ef — —^■) case, with U — 5(6) for the pe- 
riodic (disorder) case. We investigate one-dimensional 
rings of A^ = 4 and 6. Most of the results presented are 
for the A^ = 6 case, except for the results for T > where 
we have used A^ = 4 sites. The exact diagonalization cal- 
culations employ periodic boundary conditions. 

A. Ordered Clusters 

We have investigated the ground-state properties as 
a function of the hybridization and the energy spacing 
in the conduction band, A = 4t/(A - 1) = ^. We 
have calculated the average /— and c— local moments, 
< {^,{f >^< S{'^S{^^ > and < (m^^ >^< Srsr >, 
respectively. Here, S*/ and Sf are given by 



and 






^?4e- 



(2) 



(3) 



where t are the Pauli matrices. 



We have also calculated the zero-temperature /-/ and 
f-c spin correlation functions (SCF) < Sf Si_^_^ >=< 
g\Sf'S{l\\g > and < S{ Sf >=< g\S{^'S^^'\g >, respec- 

tively. Here, \g > is the many-body ground state and S'j 
is the z-component of the f-spin at site i. As expected, the 
cluster has a singlet ground state (Sg — where Sg is the 
ground-state spin) at half filling. We compare the onsite 
Kondo correlation function < Sf Sf > and the nearest- 
neighbor RKKY correlation function < Sf SJ^j^^ > to as- 
sign a state to the Kondo or RKKY regimes, in analogy 
with mean field treatments^''. The spin structure factor 
related to the equal-time f — f spin correlation functions 
in the ground state is 






S^^ |g > e»'?(-- 



(4) 



The temperature-dependent local f-spin susceptibility, 



e^ — —U^/2— -3 (symmetric case) and a mixed-valent 
(MV) type B atom with e^ = -2 and U^ = I. For both 
types of atoms V"^ ^ V'^ ^ V ^ 0.25. For t = 1, this 
choice of parameters leads to a degeneracy between the 
doubly-degenerate c-energy levels, tk = —t, and the en- 
ergy level e? -t- U^ . Upon filling the single particle energy 
levels for any x, N ~x {x) electrons fill the e4 (e?) levels, 
and two electrons the -2i conduction energy level, with 
the remaining N — 2 electrons accommodated in the x+A 
degenerate states at —t. This in turn results in strong 
charge fluctuations. 

The temperature-dependent / susceptibility, xl (T) , at 
concentration x, is given by 






1 

Q 



C^,ac^ 



<acj^^«5^°*|ac. >,(7) 



where Cx denote the configurations at concentration x, 
\ac-c > and Eq^ are the configuration-dependent exact 
many-body eigenstates and eigenvalues, respectively, and 
Q denotes the partition function. 



where 



Q = E' 



(6) 



III. RESULTS AND DISCUSSION 



Ordered Clusters 



is the partition function. Here, 5'-^°* is the z-projection of 
the total spin (both the /- and c-contributions), and \a > 
and Ea are the exact many-body eigenstates and eigen- 
values, respectively. When V = 0, the locahzed spins and 
conduction electrons are decoupled and X (T) is simply 
the sum of the Curie term due to the free f spins and the 
Pauli term of the free conduction electrons. For finite V, 
X^ {T) decreases with temperature at low-temperatures. 
The specific heat is calculated from the second derivative 
of the free energy F, Cy = — T§^. At F = 0, the specific 
heat of the system is given by the sum of the delta func- 
tion at T = that originates from the free localized spins 
and the specific heat of free conduction electrons. For 
finite V the specific heat exhibits a double-peak struc- 
ture: the high-temperature peak is almost independent 
of the hybridization and arises from the free conduction 
electron contribution, whereas the low-temperature peak 
varies strongly with hybridization. 



1. Ground State Properties 

In Fig. 1 we present the variation of the local Kondo 
SCF < S'/S'f > (squares) and the nearest-neighbor 
RKKY SCF < S(S(^^ > (circles) as a function of hy- 
bridization for two values of the hopping matrix element 
t — 0.2 (closed symbols) and t = 1.2 (open symbols), 
respectively. As expected, for weak hybridization V 
the local nearest-neighbor RKKY (Kondo) SCF is large 
(small), indicating strong short-range antiferromagnetic 
coupling between the / — / local moments, which leads to 
long range magnetic ordering for extended systems. As V 
increases, < Sl Si_^^ > decreases whereas the < Sf Sf > 
increases (in absolute value) saturating at large values of 
V. This gives rise to the condensation of independent lo- 
cal Kondo singlets at low temperatures, i.e., a disordered 
spin liquid phase. For large V the physics are local. Inter- 
estingly, as t or A decreases the f-c spin correlation func- 
tion is dramatically enhanced while the /-/ correlation 
function becomes weaker, indicating a transition from 
the RKKY to the Kondo regime. In Fig. 2 we present 
the average local /- (circles) and c- (squares) moments 
as a function of hybridization for two values of the hop- 
ping matrix element i = 0.2 (closed symbols) and t — 1.2 
(open symbols), respectively. In the weak hybridization 
limit, the large on-site Coulomb repulsion reduces the 
double occupancy of the f level and a well-defined local 
f moment is formed (/i?) — 1.0 while (/i^) = 0.5. With 
Ui {— U^ or U^), to model a Kondo-type A atom with increasing V both charge- and spin- fluctuations become 



B. Disordered Clusters 

We consider a random binary alloy cluster, A^-xBx, of 
N=6 sites and different number of B atoms, x = 0-N, ar- 
ranged in a ring described by the half-filled [N^i = 2N = 
12) two-band lattice Anderson Hamiltonian in Eq.(l). 
We introduce binary disorder in the /-orbital energy e\ 
(= e4 or e?) and in the intra-atomic Coulomb energy 



0.0- 
-0.1- 

m 

C -02- 


\ k ° p 




o 




ra -0-3- 


\r V 




0) 


yt=0.2 h t=1.2 




^ -0.4- 

o 

O -05- 
C 
■q. -06- 




<SfiS'^i> 


(n 


\ ■'■'■"■ 




C -0.7- 


■ ■'■ 




Q. 


■■'■' 




U) -08- 







0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



FIG. 1: Nearest neighbor f-f spin-spin correlations (circles) 
and on-site f-c spin-spin correlations (squares) as a function of 
V for two values of the hopping parameter of f = 0.2 (closed 
symbols) and t — 1.2 (open symbols), respectively. 



FIG. 2: /- (circles) and c— (squares) local moment versus 
hybridization for two values of the hopping parameter of t=0.2 
(closed symbols) and t=1.2 (open symbols), respectively. 



enhanced and the local /— moment decreases monotoni- 
cally whereas the c— local moment exhibits a maximum. 
In the large V limit both the /— and c— local moments 
show similar dependence on V, with < /i^ >«< /i? >, 
indicating that the total local moment /i vanishes. The 
effect of lowering the energy spacing A is to decrease (in- 
crease) the /— (c— ) local moment, thus tuning the mag- 
netic behavior of the system. Note that the maximum 
value of the c— local moment increases as A decreases. 
This is due to the fact that for smaller t values the kinetic 
energy of conduction electrons is lowered, allowing con- 
duction electrons to be captured by f electrons to screen 
the local f moment, thus leading to an enhancement of 
the local c— moment. In Fig. 3 we present the en- 
ergy spacing versus V zero-temperature phase diagram of 
the nanocluster, which illustrates the interplay between 
Kondo and RKKY interactions. In the RKKY region 
< <S'/S'/+i > is larger than the < S{ Sf > and the total lo- 
cal moment is non zero; in the Kondo regime < S^ SL-^ > 

is smaller than the < S^ Sf >, the total local moment 
vanishes, and the ground state of the system is composed 
of independent local singlets. The solid crossover curve 
indicates the F = 14 or A = Ac values, where the lo- 
cal and non local spin correlation functions are equal, 
i.e., < S{s{_^_i >=< S{ Sf >. The dashed curve denotes 
the set of points where the on-site total local moment 
fj, = 0. Thus, in the intermediate regime, which will be 
referred to as the free spins regimeii, < S( S(_^i > is 

smaller than the < S'/S'f >, the / moment is partially 
quenched and /J. 7^ 0. Interestingly, we find that the free 
spins regime becomes narrower as the average level spac- 
ing A is reduced. This result may be interpreted as a 
quantum critical regime (QCP) for the nanoring due to 



the finite energy spacing, which eventually reduces to a 
quantum critical point when A — > 0. Fig. 4 shows the 
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FIG. 3: Energy spacing A versus hybridization zero- 
temperature phase diagram. The solid curve denotes the 
crossover point of the spin-spin correlation function in Fig.l; 
the dashed curve denotes the set of points where the on-site 
total moment square ((/x/ -I- /ic)^) = 0.0 ± 0.05. 

spin structure factor of the local f electrons S^^{q) for 
various values of V and for t = 0.2. As discussed earlier, 
the ground state of the half- filled symmetric periodic An- 
derson model is a singlet. For small V, the spin struc- 
ture factor exhibits a maximum at g = tt, indicating the 
presence of strong antiferromagnetic correlations between 
the local f moments, consistent with the large values of 



< Sf Si_^^ > in Fig. 1. With increasing hybridization, 
the maximum of S^^ {q = tt) decreases and vanishes at 
very large hybridization, indicating that the ground state 
undergoes a transition from the antiferromagnetic to the 
nonmagnetic Fermi hquid phase. This is consistent with 
the zero-temperature phase diagram in Fig. 3. The spin 
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FIG. 4: Spin structure factor as a function of wave-vector for 
different values of V and for t=0.2. 



gap as a function of hybridization V for two values of en- 
ergy spacing is shown in Fig. 5. The spin gap is defined 
as the energy difference between the singlet ground state 
and the lowest-lying excited triplet (5 = 1) state. As 
expected, there is a nonzero spin gap for the half-filled 
Anderson lattice model, which increases with hybridiza- 
tion. Interestingly, the spin gap dramatically increases 
as the average energy level spacing A is reduced. Thus, 
the energy spacing or equivalently the size of the cluster 
tunes the low-energy excitation energy which controls the 
low-temperature specific heat and susceptibility. 




FIG. 5: Spin gap as a function of V for t — 0.2 and 1.2. The 
spin gap increases exponentially (linearly) for small (large) V. 



2. Thermal Properties 

The T=0 exact diagonalization results on small 
clusters are generally plagued by strong finite size 
effect32SiS£. Performing calculations at T > gives 
not only the thermodynamic properties of the system, 
but most importantly diminishes finite-size effects for 
(fcBT»A). 

In Fig. 6, we show the nearest-neighbor f-f spin-spin 
correlations and on-site f-c spin-spin correlation as a func- 
tion of temperature for for t — 0.2 and for V = 0.2 < Vc 
and V = 0.4 > Vc, where Vc — 0.25. At high tempera- 
tures, the free moments of the f and conduction electrons 
are essentially decoupled. The nearest-neighbor non local 
spin correlation function falls more rapidly with T than 
the on-site local f — c spin-spin correlations, indicating 
that the non local spin correlations can be destroyed eas- 
ier by thermal fluctuations. For V < Vc, the nanoclus- 
ter is dominated by RKKY (Kondo) interactions at tem- 
peratures lower (higher) than the crossover temperature, 
"^RKKY' which denotes the temperature where the non 
local and local interaction become equal in the nanoclus- 
ter. In the infinite system this temperature would denote 
the ordering Neel temperature. On the other hand, for 
V > Vc the RKKY and Kondo spin correlation functions 
do not intersect at any T, and the physics become domi- 
nated by the local interactions. In Fig. 7 we present the 
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FIG. 6: Nearest-neighbor f-f and on-site f-c spin-spin cor- 
relation functions versus temperature for t = 0.2 and for 
y = 0.2 < Vc and V = 0.4 > Vc, Vc = 0.25. 



crossover temperature T^kky ^'^^ t;he cluster as a func- 
tion of hybridization for different values of t. This repre- 
sents the phase diagram of the strongly correlated nan- 
ocluster, which is similar to the "Doniach phase diagram" 
for the infinite Kondo necklace model. The phase within 



the crossover curve denotes the regime where the non lo- 
cal short-range magnetic correlations are dominant. For 
V < Vc and T >> Tf^j^j^y o^^^ enters into the disor- 
dered "free" local moment regime. On the other hand, 
for V > Vc and at low T, the nanocluster can be viewed 
as a condensate of singlets, typical of the Kondo spin- 
liquid regime. Interestingly, the Tjjjf ifv '^^^ ^^ tuned by 
the energy spacing A or the size of the cluster. Thus, 
increasing A or decreasing the size of the nanocluster re- 
sults to enhancement of the non local nearest-neighbor 
magnetic correlations and hence T^^j^y- This result is 
the first exact "Doniach phase diagram" for a nanoclus- 
ter. 
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FIG. 7: Effect of energy spacing, A = j^^ on the exact 
"Doniach phase diagram" for a strongly correlated electron 
nanochain. The crossover curve represents the crossover tem- 
perature TJjif^y, where the non local short range AF spin 
correlations become equal to the local on-site Kondo spin cor- 
relations. 

In bulk Kondo insulators and heavy- fermion systems, 
the low-T susceptibility and specific heat behavior is de- 
termined by the spin gap, which for the half-filled An- 
derson lattice model, is determined by the ratio of V to 
U. On the other hand, strongly correlated nanoclusters 
are inherently associated with a new low-energy cutoff, 
namely the energy spacing A of the conduction electrons. 
Thus, a key question is how can the low-temperature 
physics be tuned by the interplay of the spin gap and the 
energy spacing. In Fig. 8 we present the local f magnetic 
susceptibility as a function of temperature for t = 0.2 and 
for T/ = 0.2 < Vc, V = Vc = 0.25, and V = 0A> Vc- For 
small V, the spin gap which is smaller than A controls 
the exponential activation behavior of x^ at low T. On 
the other hand, in the large V limit, the spin gap becomes 
larger than A (see Fig. 5) and the low-T behavior of the 
susceptibility shows no exponential activation. At high 
T we can see an asymptotic Curie- Weiss regime, typical 
of localized decoupled moments. In Fig. 9, we present 
the specific heat as a function of temperature for V = 0.4 



FIG. 8: Local f magnetic susceptibility as a function of tem- 
perature for t = 0.2 and with V = 0.2 < Vc, V = V^ = 0.25 
and V" = 0.4 > Vc. 



and different t. At V = 0, the specific heat is given by 
the sum of a delta function at T = for the localized 
spins and the specific heat of free fermions. As expected, 
by switching on the coupling V, they are combined to 
form a two-peak structure. The broad peak at high T 
is rather similar to the free-electron gas. The low-T be- 
havior is associated with the lowest energy scale, which 
as in the case of the susceptibility, is determined by the 
lowest value between the spin gap and the energy spacing 
A. For large values of t (or A) the spin gap is reduced 
(see Fig. 5) and the spin gap is the lowest energy scale. 
Consequently, the low-T behavior exhibits exponential 
activation associated with the spin gap. On the other 
hand, for small energy spacing the physics become local 
(Kondo regime) and the low-T sharp peak shifts towards 
higher temperatures and becomes broader. 



B. Disordered Clusters 

1. Effect of Disorder 

The configurations for a; < 3 are shown in Fig. 10, left 
panel, along with the value of the spin, Sg, of the ground- 
state. The A (B) atoms are denoted by closed (open) cir- 
cles, respectively. Except for the homogenous cases {x—0 
and a; = 6), with a Sg = ground state, for all x there 
are configurations with Sg ^ 0. The average occupation 
and average LM for the periodic Kondo and MV lattices 



are < n 

,B\2 



A\2 



/>=!'< W) 



> 



0.99, and < nf > 



1.6, 



< (a*/ ) > = 0.43, respectively. We carry out a detailed 
analysis for x=l (Sg =2) to demonstrate the FM transi- 
tion induced by a single MV atom in an otherwise Kondo 
cluster. Studies of extended systems have reported simi- 
lar occurrence of ferromagnetism in the MV phaseSS. As 
expected, the singlet ground state of the cc = Kondo 









■ ■ ■ ■ t=0.2 . • . 








1=0.6 , ■ / 


1.0x10"'- 


' A 




v'-^.'° ■••■■■ / 


c„ 


■A 




/ 


5.0x10"'- 


1 1 

. \ 
1 \ 


\ 


v_^^ 


0.0- 


!/ 


^ 


1 ' 1 ' 1 ' 



T(K) 



FIG. 9: Specific heat as a function of temperature with V = 
0.4 and various values of t = 0.2, 0.6 and 1.0. The low-T peak 
for larger energy spacing is due to the spin gap. 



cluster is characterized by n.n. anti-ferromagnetic (AF) 
/-/spin correlations (< Sf{i)Sf{i -|- 1) > = - 0.58). The 
introduction of a MV atom renders them ferromagnetic. 
Since Ub is small, the B impurity tends to remove charge 
from the the conduction band, in particular from the k- 
state with ek = —t, which has large amplitude at the B 
site and at the opposite A site across the ring. Such a de- 
pletion is different for the two spin states, thus yielding a 
maximum value for the f-moment of the MV atom. The 
/-/ spin correlation function between the Kondo and MV 
atoms are AF (< Sf[i)Sf{i + !)>=- 0.23), while they 
are FM among the Kondo atoms (< Sf{i)Sf{i + l) > 
= +0.94). A similar result was recently found in ab ini- 
tio calculations—, where introducing a nitrogen impu- 
rity in small (1-5 atoms) Mn clusters induces ferromag- 
netism via AF coupling between the N to the Mn atoms, 
whilst Mn-Mn couple ferromagnetically. We find that 
there is a crossover in Sg from ^ 1 ^ 2 ^ (Fig. 10, 
right panel) indicating a reentrant nonmagnetic transi- 
tion around es = 2. This almost saturated FM Sg — 2 
domain is robust against small changes in Ub, V , ^a, Ua, 
cluster size (A'' = 7), and band filling (A^e; = 10) provided 
that the Kondo atom has a large LM. 

In Fig. 11 we present Txl{T) as a function of temper- 
ature for different x. As T ^ (inset Fig. 11) Txl{T) 
approaches a finite value for x = I — A while it vanishes 
exponentially for a::=0, 5 and 6. This is due to the fact 
that the former concentrations involve some configura- 
tions which are magnetic, while the latter have singlet 
ground states (Fig. 10). The stronger (weaker) low- 
temperature dependence ior x — I {x — 2 — 4) is due to 
the smaller (larger) spin gap between the ground state 
and the lowest excited states. The magnetic suscepti- 
bility displays also a magnetic crossover upon varying x, 
and reveals a Curie- like divergence at low T for x = 1 — A. 
The temperature-dependent results for the specific heat. 
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FIG. 10; Left panel: Alloy configurations for various concen- 
trations a; < 3 (the a; > 3 cases are obtained by exchanging 
closed and open circles). For each x < 3 configuration, the 
value of the ground-state spin Sg is reported. Right panel: 
Energy difference (in units of 10~*t ) between the lowest S < 2 
eigenstates and the ground state as function of es . 
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FIG. 11: Temperature dependence of the average f- 
susceptibility for different alloy concentrations. The inset 
shows the low-temperature behavior. 



not reported here, show corroborative evidence of this 
disorder-induced magnetic crossover. 

2. Effect of Energy Spacing 

Next we address a number of important open issues, 
namely (1) the effect of A on the interplay between 
RKKY and Kondo interactions in disordered clusters, (2) 
the characterization of the single-impurity "Kondo cor- 
relation energy" Tk in a dense-impurity cluster and (3) 
the effect of disorder and A on the distribution of the 
local Tx's. In the following, e^ = —2. 

In contrast with previous studies, which introduced a 
phenomenological distribution P(Tk) of single-impurity 
Kondo temperatures, the advantage of the present cal- 
culations is that one calculates exactly the Kondo cor- 
relation energy: we employ the so-called "hybridization" 



approachai, with Tk defined as 

kBTK{l)=Eg{V,^{})-Eg 



(8) 



where Eg[Vi — Q) is the ground-state energy of the dense- 
impurity cluster when V is set to zero at the \th site. 
Eq.(8) reduces to UbTk = ^6Q„d-£^F+e/-£^;r2i2^in the 
single impurity case. Here, Ep is the highest occupied 
energy level in the conduction band and Ei,and is the 
conduction band energy. This definition of the local Tk 
takes into account the interaction of the /-moment at 
site i with the other /-moments in the system'^^. 

In Table I we list for the periodic, x—Q, case the local 
Kondo /-c spin correlation function < Sf{i)S^{i) >, the 
n.n. /-/spin correlation function < Sf{i)Sf{i + 1) >, 
and the local Kondo temperature for two different values 
of t (The energy spacing is A = 4t/(7V - 1) = 4t/5). 
As i or A decreases the f-c spin correlation function is 
dramatically enhanced while the /-/ correlation function 
becomes weaker, indicating a transition from the RKKY 
to the Kondo regime. This is also corroborated by the 
increase in the local T^(j). The energy spacing affects 
not only the magnetic (A) atoms but the MV atoms as 
well. Thus, increasing t drives the B atoms from the 
non-magnetic, NM {nj « 2), to the MV and finally to 
the Kondo regime. 

We next examine the role of even versus odd num- 
ber of electrons on the magnetic behavior of the uni- 
form x=0 case. For t = 1, changing the number of elec- 
trons from Nej = 12 to N^i = 11 results in: (a) an en- 
hancement of the local Kondo f-c spin correlation func- 
tion, < Sf{i)S^ii) > from -0.01 to -0.12; and (b) a 
suppression of the nearest-neighbor f-f spin correlation 
function < Sf{i)Sf{i + I) > from -0.58 to -0.20 (Due 
to the broken symmetry for Ngj = 11, the f-f spin cor- 
relation functions range from -0.5 to -1-0.02). This in- 
teresting novel tuning of the magnetic behavior can be 
understood as follows: For an odd-electron cluster, the 
topmost occupied single particle energy level is singly 
occupied. On the other hand, for an even-electron clus- 
ter, the topmost occupied single-particle energy level is 
doubly occupied, thus blocking energy-lowering spin-flip 
transitions. This energy penalty intrinsically weakens the 
Kondo correlational^. As expected, changing the number 
of electrons from even to odd changes Sg = to S'^ = i . 
For t= 0.05 (Kondo regime), the on site f-c correlation 
function does not depend as strongly on the parity in the 
number of electrons because the sites are locked into sin- 
glets. On the other hand, < Sf{i)Sf{i 4- 1) > becomes 
suppressed as in the case of large energy spacing. Similar 
results were found for the various disordered concentra- 
tions. 

In Fig. 12 we present the local Tff(i) as a function of 
the local f-c spin correlation function < Sf{i)S^{i) > 
for all Kondo (A) atoms in the singlet ground state at 
any concentration x for t= 0.05 and 1.0. Note the differ- 
ent scales both on the horizontal and vertical axis in the 
panels. In both panels, the closed circles correspond to 
the x=0 lattice case and the line is a guide to the eye. 
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FIG. 12: A-atoms: Local Kondo Temperatures (in K) vs the 
local f-c spin correlation function, for different configurations 
and two different values of t. The closed circles refer to the 
X — case and the lines are a guide to the eye. 



TABLE I: Local Kondo f-c and n.n. f-f spin correlations 
functions and the local Kondo temperature (in K) for two 
values off (in eV). The average energy spacing is A = 4t/{N — 
1) = 4f/5. 





< Sf{i)S^ii 


) > <Sf{i)Sf{i + l) > TK{i) 


t=0.05 


-0.626 


-0.322 551.8 


t=1.00 


-0.011 


-0.584 173.4 



The results indicate a correlation between Tk and the 
f-c spin correlation function (the larger T/f 's correspond 
to the more negative f-c values) as one would expect, 
since both provide a measure of the Kondo effect. For 
t=0.05, most of the disordered cluster configurations are 
in the Kondo regime {Sg = 0), with larger Tk values; 
consequently, panel (a) has a larger number of singlet 
configurations. The introduction of MV impurities in- 
duces a distribution of Tr-(j)'s, whose values are overall 
enhanced compared to those for the x=0 case, except for 
several configurations for i=0.05, in contrast with single- 
site theories for extended systems^. It is interesting that 
P{Tk) for t=0.05 exhibits a bimodal behavior centered 
about 710 and 570K, respectively: The higher T^'s origi- 
nate from isolated Kondo atoms which have MV atoms as 
n.n. so that the local screening of the magnetic moment 
of the A atom is enhanced. 

The effect of alloying and A on the RKKY versus 
Kondo competition for a given x is seen in Fig. 13 
(left panel), where the configuration averaged local < 
Sf{i)Sf{i) >x and < Sf{i)Sf{i -\- 1) >j: correlation 
functions arc plotted as a function of t. The solid curves 
denote the uniform x—0 case, where we drive the clus- 
ter from the RKKY to the Kondo regime as we decrease 
t. We find that the stronger the average Kondo correla- 
tions are the weaker the average RKKY interactions and 
vice versa. In the weak Kondo regime the configurations 
exhibit a wider distribution of RKKY interactions indi- 
cating that they are sensitive to the local environment. 
In contrast, in the strong Kondo regime, the Kondo (A) 
atoms become locked into local Kondo singlets and the 
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FIG. 13: Left panel: Configuration-averaged local f-c (top) 
and n.n. f-f spin correlations (bottom) for the A atoms as 
function of f . The solid line refers to the homogenous a; = 
case. Right panel: Zero-temperature t vs a: phase diagram 
for the nanocluster. Black (gray) circles denote the Kondo 
(RKKY) regime. The white circles and the dashed contour 
delimit the FM region. The horizontal stripes denote the non- 
magnetic (NM), mixed valence (MV) and local moment (LM) 
behavior of the B-atoms. 

Fig. 13 we present the t versus x phase diagram for the 
nanocluster at T = . We compare the < Sf{i)S^{i) >x 
and < Sf{i)Sf{i + l) >a; to assign a state of specific con- 
centration to the Kondo or RKKY regimes (black and 
gray circles, respectively), in analogy with the a; = case 
(Table I) and with mean field treatmentsSi. The hori- 
zontal gray stripes denote qualitatively ranges of t where 
the B atoms exhibit NM, MV and LM behavior. An in- 
teresting feature of the phase diagram is the appearance 
of a large FM region {Sg ^ 0) enclosed by the dashed 
line. The RKKY region at large t and large x originates 
from the B atoms which become magnetic. For the non 
FM configurations and for x < 5 the Kondo (RKKY) 
correlations of the A atoms dominate at small (large) t, 
in analogy with the a; = case. On the other hand, for 
X = 5 the local Kondo correlations of the single A atom 
at low t dominate over the /-/ correlations between the 
A-B and B-B pairs. For the uniform {x=6) MV case 
we include only results in the large t regime, where the 
MV atoms acquire LM's which couple antiferromagnet- 
ically. Overall, the RKKY interactions prevail for any 
concentration when t is comparable or larger than the 
hybridization V. 



IV. CONCLUSIONS 

Recent advances in STM experiments have made it 
possible to study the electronic and magnetic properties 



of strongly correlated electrons in nanoscopic and meso- 
scopic systems. There are two main differences between 
nanosized clusters and the infinite lattice. First, the dis- 
crete energy levels of the conduction band states intro- 
duce a new low-energy scale, i.e., the average energy level 
spacing A. This new energy scale that competes with the 
spin gap can effect the low-temperature behavior of the 
system. Second, the results depend on the parity of the 
total number of electrons. li Ntot is odd, the ground state 
is doubly degenerate. 

We have carried out exact diagonalization calculations 
which reveal that the: (1) energy spacing; (2) parity of 
the number of electrons; and (3) disorder, give rise to a 
novel tuning of the magnetic behavior of a dense Kondo 
nanocluster. This interesting and important tuning can 
drive the nanocluster from the Kondo to the RKKY 
regime, i.e. a tunable "Doniach" phase diagram in nan- 
oclusters. We have employed the criterion of comparing 
the exact non local versus local spin correlation functions 
to determine if the nanocluster lies in the RKKY ver- 
sus Kondo regime. For weak hybridization, where the 
spin gap is smaller than A, both the low-temperature 
local / susceptibility and specific heat exhibit an ex- 
ponential activation behavior associated with the spin 
gap. In contrast in the large hybridization limit, A is 
smaller than the spin gap, the physics become local and 
the exponential activation behavior disappears. The in- 
terplay of A and disorder produces a rich structure zero- 
temperature alloy phase diagram, where regions with pre- 
vailing Kondo or RKKY correlations alternate with do- 
mains of FM order. The distribution of local Tk and 
RKKY interactions depends strongly on the local envi- 
ronment and are overall enhanced by disorder, in con- 
trast to the hypothesis of single-impurity based "Kondo 
disorder" models for extended systems. We believe that 
the conclusions of our calculations should be relevant to 
experimental realizations of small clusters and quantum 
dots. For example, the recent experiments^^ of magnetic 
clusters on single- walled carbon nanotubes of varying size 
provide much flexibility for investigating the interplay of 
Kondo and RKKY effects at different energy scales. 
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